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Abstract 
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We study the evolution of nearest-neighbor entanglement in the one dimensional 
Ising model with an external transverse field. The system is initialized as the so 
called "thermal ground state" of the pure Ising model. We analyze properties of 
generation of entanglement for different regions of external transverse fields. We find 
that the derivation of the time at which the entanglement reaches its first maximum 
^-j- I with respect to the reciprocal transverse field has a minimum at the critical point. 

OO ■ This is a new indicator of quantum phase transition. 
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1 Introduction 

Entanglement has been recognized as an important resource for quantum information and 
computation pQ, and has been the subject of intense research over the past few years [2]. 
Recently, there has been extensive analysis of entanglement in quantum spin models [3]. 
Various models have been considered for entanglement generation, and their static [4, 5J 
as well as dynamical properties [6j [7J |8] have been investigated. The reason for this is 
that, many of these models can be realized in cold atomic gas in an optical lattice [9| [TO]. 

A special one among these models is the anisotropic XY model. It can be solved 
exactly by means of the Jordan- Wigner transformation [llj. A lot of works have been 
done on the dynamics of entanglement in this model. In Ref. [6], the time evolution of 
initial Bell states was studied. Ref. [7] investigated the dynamics of entanglement of a 
system prepared in a thermal equilibrium state, and the system starts evolving after the 
transverse field was turned off. 
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In this paper, we analyze a spin- 1/2 transverse Ising chain with a particular separable 
initial state, the so-called "thermal ground state" pi] for the case of zero field. We study 
the generation of nearest-neighbor entanglement as a function of the external transverse 
field, and find the time at which the entanglement reaches its first maximum. It turns 
out that the derivation of this time with respect to the reciprocal transverse field has a 
minimum at the critical point. This is a new indicator of quantum phase transition. 



2 Static correlations 

The anisotropic XY model is described by the Hamiltonian 

\ N i N 

H = "i Ek 1 + ^ a ?+i + 7K°i+i] - 2 E °? ' w 

i=\ i=l 

where <rf are the Pauli matrix of the spin at site i and we assume periodic boundary 
conditions. We have neglected a common factor and absorbed the external field into the 
reciprocal field A. The anisotropy parameter 7 connects the quantum Ising model for 
7 = 1. In the range < 7 < 1, the model belongs to the Ising universality class and 
undergoes a quantum phase transition at the critical point A c = 1. We will consider the 
case of 7 = 1 in this paper, namely the transverse Ising model 

1 - 

#=-^E( A «+i + ^) • (2) 

i=l 

The Hamiltonian defined in Eq.(2) can be mapped into a one-dimensional spinless 
fermion system with creation and annihilation operators cj and q via the Jordan- Wigner 

transformation '"^ = Yij<i( l ~ 2c j c i) c i 5 ~ y ^ = 4 Tij<i( l - 2c i c i)' \ = 4 c i ~~ §■ By 
making use of the transformation 

^ k = -jfi E ^("^ + iPk4) ( 3 ) 
further, with ot^ = A fc -(i+Acosfc) n _ Asmfc we ^ ^ e fermionic 

^^-(l+Acos^Afe] ' ^ v /2[A|-(l+Acosfc)A fc ] ' to 

Hamiltonian finally 



where the spectrum is A fc = yl + A 2 + 2A cos k. 

The evolution of the system is governed by Eq.(2). In the Heisenberg picture, the time 
evolution of the spinless fermion operator Cj(t) is [6] 

Cj(t) = — J^[cos - i)A(k, t)ci + sin fc(Z - i)B(k, t)c\] , (5) 
k,i 
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where A(k,t) = e iAkt - 2i/3| sin A fc t, B(k,t) = -2ia k f3 k sin A k t. 

The dynamics of the system also depends on the initial state. We choose the "thermal 
ground state" po = + °f the pure Ising model (no transverse field 

is applied) as the initial state. Here \N + ) = | — — >)n, and \N~) = | <— <— )n 
are the two degenerate ground states of the Ising model with all spins pointing to the 
positive (negative) x direction. Note that |iV + ) and \N~) are the ground states of the 
Hamiltonian (2) with A — > oo. When A is finite, they are even not eigenstates of the 
Hamiltonian. Thus, the evolution from this initial state must be nontrivial for finite A. 
Physically, this can be viewed as an instantaneous, i.e., idealized sudden quench [8] in the 
reciprocal field Ai — > A 2 with Ai — > 00 and A2 = A. 

Due to the translational symmetry, we need only to consider averages of the form 
(cx...)q. Namely, we choose site-1 as "the first site", where (...) denote the average over 
the initial state po- We rewrite {N^) in the fermionic representation, 

/ 1 \N 

\N±)= (-^j (l±c t 1 )(l±c 2 )...(l±4)|0 1 ,0 2 ,...,0 w ) , (6) 

where |0i, O2, ...,0^) denotes the vacuum of the fermions, or the state for all spins down. 
Here we have to pay attention to the order of operators in this equation. We write ( 1 ± c{ ) 
before operators on all other sites in order to be compatible with the correlation (ci...) . 
In other words, if another site i is chosen as "the first site", we should write the state 
as (A^) = (-i=) Ar (l ± c|)(l ± c| +1 )...(l ± 4_i)|0i, i+ i, .... 0j_i) to calculate averages like 
(cj...)o, due to the periodic boundary conditions. 

The only single-site averages are (ci)o and (4 c i)o (or their complex conjugate). It is 
easy to see that (N + \ci\N + ) = | = — (N~ |ci| N~), so (ci)o = 0. As to (c\ci)q, we have 
(N + \c\ Cl \N + ) = i(0i|(l + c 1 )c t 1 c 1 (l + 4)101) = \. Similarly, (N-\c\ Cl \N-} = §. So that, 

(4 c i)o — |- We can a ^ so see ^ ms P orn t from = (^ L )o = (4 c i — |)o- 

Next we consider two-point correlations of the form (ciq) , / > 2. For / = 2, we have 

(7V + | Cl c 2 |iV + ) = f^j (0 1 ,0 2 |(l + c 2 )(l + c 1 )c 1 c 2 (l + c t 1 )(l + 4)|0 X ,0 2 ) 

-J (0x1(1 + Cl ) Cl (l - 4)10:) • (0 2 |(1 + c 2 )c 2 (l + 4)|0 2 ) 
1 

Similarly (N~\ciC 2 \N~) = — |. So that (cic 2 ) = —\. We can show that (4 c 2)o = \ in 
the same manner. For / > 2, when q or c] crosses (1 ± 4) to its right side, this will make 
it be (1 =F 4)- But (0 2 |(1 ± c 2 )(l =F 4)1^2) = 0. So we have obtained all the two-point 
correlations 

1 



4 c i)o = t(^+1J + <*i-l,j + 2 $ 



(ciCj) = -(-Si+u + Si-tj), 

(4 c i)° = 4(^+1^ ~ h-ij). (7) 
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By the same means one can check that the averages involving three fermion operators 
all vanish. 



3 Reduced density matrix and evolution of concur- 
rence 

Making use of the periodic boundary conditions, we can focus on the reduced density 
matrix p 12 (t) at time t of the first and second spins. In the {| ft)) I tl)j I It}) I II)} basis, 
the matrix elements reads, 

p\i(t) = (n\ P 12 (t)\n) 

= (tt \p 12 (t)sts^s+s 2 \ tt> + (tl \p 12 (t)sis^s+s 2 \ ti) 
+Ut \ P 12 (t)s+s^s+s 2 \ it) + (II \ P 12 (t)s+s^s+s 2 \ u> 

= tr 12 (p 12 (t)S+S^S+S 2 ) 
= tr(p(t)SfSrS£S 2 -) 
— {^\Cic\c 2 )t- 

Other matrix elements can be obtained similarly. Thus, we have 

/ (4ci4c 2 )i -(clcic 2 ) i (4c 2 ci) t (c 2 Ci)i 



p u (t) 



(8) 



-{c\c{ci) t {c\cic 2 c\) t -(cic 2 ) t (cic 2 4) t 

(44c2)t "(c 2 4)t (Cl44 c 2)i (Cl4 c 2)t 

V (44)* (c244)< (4 c i4)< (ci4 c 24)t / 

All the elements can be equally evaluated in the Heisenberg picture, e.g., p\f(t) = 
(4ci4 c 2)t = (c\(t)ci(t)cl(t)c2(t)) . We have mentioned that all the static averages in- 
volving three fermion operators equal zero. From Eq.(5), the fermion operators in the 
Heisenberg picture are merely some linear superpositions of the same set of operators in 
the Schrodinger picture. So all matrix elements involving three operators in Eq.(8) also 
vanish. Using Wick's theorem, one gets all non- vanishing elements of p 12 (t) as following 

Pn(t) = (4(t)ci(t)) (4(t)c 2 (t)) - (c\(t)cl(t)) { Cl {t)c 2 (t)) - (4(t)c 2 (t))(4(t) Cl (t)) , 



pim = pim = (ci(t)c 1 (t)) -p\ 2 1 (t), 

Pl 2 (t) = l-p} 2 (t)-2p 2 2 (t), 
pll(t)=p 12 *(t) = (c 2 (t)c l (t)} , 



p&(t)=pg(t) = -Mt)<i(t))o. 
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(9) 



We see that if we can calculate (c\(t)ci(t))o, (c l (t)c2(t))o and (ci(t)c2(t)) , then the re- 
duced density matrix is totally determined. 

Using Eq.(5) and (7), we find, in the thermodynamic limit, that 

(c\{t) Cl {t)) = P dk sin k(A*(k)-A(k))(B(k)-B*(k)) + ^- r dk{\A{k)\ 2 + \B{k)\ 2 ) 



+ 
1 

7T 



Att 



dk cos k{\A(k)\ 2 - \B(k)\' 



dk 



A sin 2 k sin 2 A^t 
A] 



+ 



<<i(*M*)> 



- / dfc(l + cos2A;)(|A(A;)| 2 - |fi(A;)| 2 ) + / d&cos A;(|yl(£;)| 2 + \B(k)f 
dA;sin2ft[A(ft) — il*(fc)]-B*(fc) 



+ 8tt 
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A cos k sin 2 sin 2 A^i 1 



(ci(t)c 2 (t)) 



^- J dk sin2k A(k)B(k) + ^- J dk[B(k) 2 - A(k) 2 ](l - cos 2k) 

1 /' 7r 2 ism2A k t Asin 2 A fc t 

— / dk sm k[ coszAfct — 2 -= (cosk + A)J 

2vr J A fc A k 



(10) 



Note that (c\(t)c2(t)) is indeed real, so that we have p\%(t) = p}&.(t)- 

For bipartite entanglement, a commonly used measure for arbitrary states of two 
qubits is the so-called concurrence [12J. The concurrence is defined as 



C{t) = max{0,2A max (t)-trv^(tj^(t)}, 
p 12 {t) = (yy®Oyp l2 *(t)a y ®Oy, 



AD 



(12) 



where X ma x is the largest eigenvalue of the matrix a/ p 12 (t) p 12 (t) . 
In the present case, 



P 12 (t) 



( p\m 








p\l(t) \ 





pirn 


pirn 








pirn 


pim 





V P\f(t) 








pirn ) 



(13) 
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Figure 1 : The evolution of bipartite entanglement for the initial state being the "thermal ground state" 
of the Ising model. 



the four eigenvalues of a/ p 12 (t)p 12 (t) are of the form 

Ht) = \\p\l(t)\ + J Pll(t)pll(t)\, 



Ht) = ni(t)\-y/pli(t) P ll(t)\, 



Ht) = ipm+timi 



x 4 (t) = \ P \i{t) - P \i{t)\. 



(14) 



In Fig.l, we give a plot of the evolution of entanglement as a function of time t and 
parameter A. At the beginning, the reduced density matrix is 

/ 1 1 \ 



1 



P 12 (0) = 7 



(15) 



110 
110 
\ 1 1 ) 

So that the concurrence equals zero. At any later fixed instant, the entanglement as a 
function of A may have different behaviors in different time intervals. For times before 



6 




Figure 2: (a): The concurrence at different 
times t = 0.5, 1.0, 1.5, 2.0 as a function of A, 
in the region of t < 2.5, 



(b): The concurrence at different times 
t = 4.0, 5.0, 6.0, 7.0 as a function of A, 
in the region of t > 3.0. 



about t = 2.5 (Fig. 2 (a)), the entanglement maintains its magnitude for large A with 
slow oscillations. As A increases, the height of each peak deceases gradually, and finally 
vanishes as A — > oo. This is easy to understand. For A — > oo, the Hamiltonian reduces to 
the pure Ising model, and the initial state is just the eigenstate of the Ising model. Thus, 
only a phase factor contributes to the times evolving state, which does not change the 
entanglement. Physically, this corresponds to Ai = A2 = 00, hence no quench takes place. 
For the narrow band 2.47 < t < 2.69, there is no entanglement at all for all A. For times 
greater than around t = 3.0 (Fig. 2 (b)), the entanglement emerges abruptly, but only has 
non-vanishing values for very short interval of A. As the time increases, the generations of 
entanglement are different from each other for different A. For A is less than about A = 0.8 
(Fig. 3 (a)), the concurrence will encounter more than one maximums. In a time region 
between the first two maximums, no entanglement appear. Although weak oscillations 
emerge after this region, the entanglement can last for a fairly long time with a not too 
small magnitude. For 0.80 < A < 1.07, there is only one maximum as the time increases. 
After a short time, the entanglement vanishes. For A > 1 (Fig. 3 (b)), strong oscillations 
emerge. The entanglement reaches its maximum fast, then suddenly decreases to zero. 
Such a behavior will repeat again at later times, and finally the entanglement vanishes 
completely. 

For any A, there is a time T m (A) at which the entanglement reaches its first maximum 
(Fig. 4). This time is a monotonically decreasing function of A. Taking its derivative 
with respect to A, we find that the derivative dT ^ X ^ has a minimum at the critical point 
A = 1. This interesting property maybe an indicator of the quantum phase transition at 
the critical point. 
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Figure 3: (a): The dynamics of entanglement 
for different reciprocal fields A = 0.2,0.4,0.6,0.8, 
in the region of A < 0.8, 



(b): The dynamics of entanglement for 
different reciprocal fields A = f.5,3.0, 
in the region of A > f .07. 




4 Brief discussions and conclusions 



In this paper, we have studied the dynamics of the transverse Ising chain prepared in 
the "thermal ground state" of the pure Ising model. For different values of the reciprocal 
field, the evolution of the entanglement has dramatic distinctions. We also obtained the 
times at which the entanglement reaches its first maximum, which is a monotonically 
decreasing function of the parameter A. The derivative of this quantity with respect to 
the reciprocal field gets its minimum at the critical field. This is a new indicator of the 
quantum phase transition. Let us finally remark that the initial state we choose is simple 
to be prepared experimentally. One can first cool the pure Ising system down to near 
absolute zero, then turn on some desired transverse field at t — to generate the wanted 
entanglement. 
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